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ABSTRACT 



The line caustic behavior has been discussed since Chang and Refsdal (1979) 
mentioned inverse-square-root-of-the-distance dependence of the amplification of 
the images near the critical curve in a study of a single point mass under the 
influence of a constant shear due to a larger mass. A quarter century later, 
Gaudi and Petters (2001) interprets that the distance is a vertical distance to the 
caustic. It is an erroneous misinterpretation. 

We rehash Rhie and Bennett (1999) where the caustic behavior of the binary 
lenses was derived to study the feasibility of limb darkening measurements in 
caustic crossing microlensing events, (i) J = ±a/4<5u; 2 _ J_ where 5u \\ dJ, and 
<5co>2- and J_ are i?_-components of 5u (the source position shift from the caustic 
curve) and 2d J (the gradient of the Jacobian determinant) respectively; (2) The 
critical eigenvector ±E_ is normal to the caustic curve and easily determined 
from the analytic function K-field; (3) Near a cusp (J_ = 0) is of a behavior of 
the third order, and the direction of d J with respect to the caustic curve changes 
rapidly because a cusp is an accumulation point; (4) On a planetary caustic, 



\dJ\ ~ v/i/epi is large and power expansion does not necessarily converge over 

the size of the lensed star. In practice, direct numerical summation is inevitable. 

We also note that a lens equation with constant shear is intrinsically incom- 
plete and requires supplementary physical assumptions and interpretations in 
order to be a viable model for a lensing system. 

Subject headings: gravitational lensing - binary stars 



There are two normals. One is the critical eigenvector and normal to the caustic 
curve. The other is gradient J and normal to the critical curve. The angle between 
the two normals changes along the curves. When the two normals are across, the 
caustic curve stops and turns around, and the punctuation point is called a cusp. 
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l A Tale of Two Curves 
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-1. Back at Table with RB99 

The ingredients are a 2-d plane, two 2-d variables called source position and image po- 
sition defined on the 2-d plane, a relation between the source position and image position 
called the lens equation, critical curve of the lens equation where the lens equation is sta- 
tionary (or degenerate), and caustic curve onto which the critical curve is mapped under the 
lens equation. The lens equation happens to be an explicit function from an image position 
to its source position: u = f(z,z); u = f(z,z) where z is an image position and u> is its 
source position. The lens equation can be considered a mapping from the complex plane 
(z-plane because it is parameterized by z) to itself (cu-plane because it is parameterized by 
a;). This complex plane (or the underlying 2-d real plane) is referred to as the lens plane 
because that is where all the tranverse objects related to lensing are defined and studied. It 
is in fact convenient to put the lens plane at the distance of the center of the mass of the 
lensing objects especially when the focus of the lensing study is on the lensing objects such 
as micro lensing planet systems or compact objects. 

The issue is the local behavior of the lens equation in the neighborhood of the critical 
points (points on the critical curve) of quasi-analytic lens equations. A quasi-analytic lens 
equation is by defintion a lens equation that is not analytic itself but whose derivatives are 
all determined by one analytic function called /t-field and its derivatives. The n-point lens 
equations including constant shear are quasi-analytic. 

The critical curve is one of the family of equipotential curves where the potential is given 
by the magnitude of k and the variable along the critical curve is the phase angle of k. The 
phase angle also determines the direction vectors ±E± (rb99.7). The criticality defines the 
critical direction ±E_ (rb99.5), and the caustic curve is everywhere normal to the critical 
direction. An arbitrary infinitesimal deviation dz from the critical curve causes uo to shift 
only in the tangent to the caustic curve (rb99.8). The whole dimension of EL direction or the 
E_ component of an arbitrary dz is projected out by the lens equation, and that is what is 
referred to be degenerate and why the Jacobian determinant of the lens equation vanishes on 
the critical curve. In other words, 1) if we restrict to the linear order, an arbitrary deviation 
duj = 5ui + E + from a caustic point has infinitely many solutions dz = dz + E + + dz-E- where 
2dz + = 5ui + and dz^ is indeterminate; 2) thus, we need to include the next nonvanishing 
order terms to see the proper behavior of the lens equation in the neighborhood of the critical 
points. The second order terms have been written out in rb99.10 and rb99.13, or in equation 
(61) below. The third order terms are shown in equation (63). 

It is worth pointing out that the orthogonal decomposition coefficients are always real 
as should be clear from equation (51). The differential structure is about the underlying 
linearly independent two dimensional space whether parameterized by two real variables 
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or by a complex variable and its complex conjugate, and thus whether one uses real basis 
vectors (implicitly 2-column vectors (1, 0) and (0, 1) defined at one point on the caustic curve 
and the corresonding point on the critical curve) or ±E± that are tangent and normal vector 
fields on the caustic curve, we inevitably deal with the same two degerees of freedom, namely 
the tangent component and normal component. One of the obvious advantages of maximally 
utilizing the quasi-analytic nature of the quasi-analytic lens equations by employing complex 
coordinate systems is that we have explicit expressions of the basis vector fields E±. For 
example, we know how they change along the caustic curve parameterized by ip half the 
phase angle of the K-field: k — \k\ expi2ip. 

E + = e~ ilp , E_ = ie~ ilp ; dE + = -E^dip , dE_ = E + d<p . (1) 

(We note that the basis vector fields are not necessarily smooth everywhere.) 

We do see a potential problem in the emphasis to have used real variables in Gaudi & 
Petters (2001) because they are putting the focus in a worng place. If it was a methodology to 
generate an undisclosed impression on the readers against our having consistently employed 
the complex plane, it perhaps is a wronful deed. If Gaudi & Petters (2001) stated it because 
it is simply true that they worked with real coordinates, the fault lies in their failure to clear 
the smoke and address the relevant issues: 

1. intrinsic understanding of the quadratic behavior of the lens equation; 

2. correct formulation of (rb99.16 and equation(23)) and its proper interpretation; 

3. inconsistent truncation prescription in the first paragraph in SEF.p.189. 

We will address the correct issues 1 in the following, of course using complex coordinates 
(z for the image space and u> for the source space) and the analytic function /t-field. As 
corollaries, the followings will be clear. 

1. The claim by Gaudi & Petters (2001) that they reproduced SEF.6.17 is free of content. 

2. The claim by Gaudi & Petters (2001) that 5 = |Ac<j_| is a "vertical distance" is 
erroneous. 

Gaudi & Petters (2001) warns the readers, "there are not many people who know that the 
distance is not the shortest distance but the vertical distance. " The distance 5 is neither the 



Not all questions grammatically correct are right questions. 
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shortest distance nor a vertical distance as is described in the abstract and illustrated in Fig. 
1. When 5 can be approximated by the shortest distance, the "vertical distance" converges 
to 5 as well. 



Regurgitation of RB99: In the case of the binary lens in rb99.Fig.l where the two masses 
are comparable and their separation is ~ 0(1) (where 1 refers to the Einstein ring radius of 
the total mass), the gradient of the Jacobian determinant |VJ| ~ 0{1). Equation rb99.13 
reads as follows when dz + is reinstated. 

5uj^ = <5u> 2 _ = ^ {J-(dzi - dz\) + U+dz+dz-) (2) 

If we assume the linear approximation for dz + as in rb99, the quadratic equation in dz- can 
be solved in terms of dz+ and 5uj_. 




VJ\ 2 dz 2 + 4^2_ 

J 2 + ~1~ [6) 

Since dz_ is real, as we have emphasized that E± decomposition components are real variables 
above and in equation (51), dz- has solutions when the inside of the square root is non- 
negative (rb99.14). 

5uo-d-J + \dJ\ 2 dz 2 + > (4) 

The degenerate solution for which the equality of equation (4) holds amounts to a shift along 
the critical curve. 

dz_ = — y-dz + =3- J + dz + + J_dz_ = (5) 

If z is a critical point, then z' Q = z Q + dz + (E + — J + E_/ J_) is also a critical point. 

We rewrite z Q + dz as a vector deviated from z' such that z Q + dz = z' Q + dz'. Then, dz' 
is a vector in the direction of ±E^(z Q ). 



. ,. , , . , \VJ\ 2 dzl 45io 2 - / _\ 

dz'(z' ) = ±E_(z ) ' J2 + + —j- (6) 

If u = u(z ) is the caustic point corresponding to the "old" critical point z , then uj' = 
u + 2dz + E + (z a ) is a new caustic point because E + (z Q ) is tangent to the caustic curve. We 
can easily calculate the new caustic point u' Q up to the second order. Let's start for clarity 
with the general second order expansion in dz in equation (48). 

5uo 2 = \dRdz 2 = -(-dJ)(dz+ - idz_) 2 (7) 
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The second equality is obtained by using 3 J = —kBr, .^-decomposition dz = dz + E + + 
dz_E_, and miscellaneous relations such as E_ = iE + that can be found in rb99 and 
section 2. In order to get the expression for the caustic point, we apply the condition 
= J + dz + + J_dz_ and get 5u parameterized by dz + . 

( 1IV/I 2 \ 1IV/I 2 

5u = 5uo + E + + SuS. = { 2dz + + - l —j^J + dz 2 + J E + - -~^J-dz 2 + E_ (8) 

We rewrite the equations for 5u + and 5uo^ as follows for later use. 

= {8u+-2dz+)J- + J+8u- (9) 



1 IVJ 



2 



5lu_ = -ILLU-dzl (10) 

4 J_ + 



It should be clear that the role of dz + is mainly to cause a migration of the cautic point 
and critical point from oj and z to new ones u)' and z' respectively. In order to see the 
uncluttered picture of the local behavior of the lens equation near a critical point, we set 
dz+ = 0. 



-1. In Details We Trust: Brutally Honest Interpretations 

If dz + = 0, dz = dz-E-, and from equation (61), 

5uj 2 + = -dz 2 _J± 
4 

5uj 2 +E + + 5uj 2 -E_ = -dJ dz 2 _ . 

This is the quadratic behavior of the lens equation near a critical point. The role of the 
critical direction is clear. If z is a critical point, two images separated along the critical 
direction at z — z Q ±\dz_\E_ are from a same source. If u is the caustic point corresponding 
to z under the lens equation, then the position of the source, uj = u + 5lj, that produces the 
two images is such that the shift Sou is in the same direction as the gradient of the Jacobian 
determinant VJ = (J+, J ) at z . See Fig. 1. We refer to this caustic point u as the 
preferred caustic point for uj. 



5u)\± = ; 



Question : Suppose uj is an arbitrary source position near the caustic curve. Then, it is 
conceptually straightforward how to find the images of uj near the critical curve. 



- 6- 



Solution : We find a caustic point uj such that 5uj = uj — uj q is parallel to dJ (orV J) at z Q 
where z Q is the critical point that corresponds to uj . Then, two images of uj can be found 
in the critical direction from the critical point z Q : z = z Q ± \dz-\E- where 



ldz - l = \l— = \l— = yw 

The two images have opposite parities and same magnification as reflected in the signs and 
absolute value of the Jacobian determinant at the image positions (rb99.11). 

J(z + dz) = d.J = dz_J_ = ±|dz_||J_| = ±^A5uj 2 -J- (12) 

If Suj is antiparallel to dJ, the inside of the square roots in equation (11) is negative, and uj 
generates no images near the critical curve. 

1. Once the directionalities of the image line (±E_) and source line (dJ) are determined, 
the lens equation becomes a real symmetric quadratic equation from the image line 
to the source line. The lens equation from the image line to the source line has two 
solutions for 5udJ > (real and positive), no solutions for 5uj8J < (real and 
negative), and degenerate at the critical point where 5ujdJ = 0. It is exactly like the 
simple real quadratic equation y = ax 2 : y = ax 2 has two solutioins for ay > 0, and 
no solutions for ay < 0. At the critical point x = where dy/dx = 0, the solutions 
are degenerate. 

It should be worth emphasizing that the image line (±E_) and source line (8 J) are 
not orthogonal to each other even though it may be tempting to imagine so because of 
our perception from the graph y = ax 2 which we are accustomed to draw in orthogonal 
grids. 

2. In fact, the second order approximation fails where the image line and the source line 
are orthogonal (J_ = 0). Since equation (3) fails where J_ = 0, we need to consider 
the third order terms. The underlying reason can be seen in that the second order 
terms accomodate only a migration along the critical curve as is clear from equation 
(12). 

3. Fig. 1 clearly shows that the larger the angle between VJ and E_, the more tangen- 
tially the images z± are shifted from z Q , and the smaller the normal distance of the 
images from the critical curve which determines the J- value (dJ). If we consider a 
small stellar disk near the caustic curve, the images will be extensively elongated along 
the critical curve where J+ > J_. 
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4. Geometrically, J_ = means that the critical direction is tangent to the critical curve. 
Let's consider the lens equation as a mapping restricted to the critical curve. Then, 
the restricted mapping is stationary or critical where J_ = 0. In other words, the 
caustic curve as an integral curve of the tangent develops a point (cusp) where J_ = 
because the "speed" vanishes as in the trajectory of a pendulum 2 at the maximum 
gravitational potential. 

dcu = duo + E + + duo„E_ = / 2dz + E + (13) 

J J J critical curve 

= dJ = dz + J + + dz-J- = dz + J + =>- 5lu + = 2dz + = (14) 

The cusps are the only singularities on the caustic curve. That is because the lens 
equation is smooth on the critical curve, and the critical curve is smooth. (A critical 
curve can have bifurcation points that are four-prong vertices, and the caustic curve 
develops cusps; see RBCPLX). Thus, the caustic curve is smooth except at the critical 
points where the restricted mapping is stationary. 

5. If J + = 0, the two normals to the critical curve and the caustic curve coincide. 
|J_| = |VJ| and |5cj 2 -| = \Su\. The image line (the line connecting the two im- 
ages) is perpendicular to the critical curve, and \5u\ is the (shortest) distance from oo 
to the caustic point onto which the critical point that disects the image line is mapped 
under the lens equation. There is no need to be an alarmist for that the second term 
in equation (11) is not well defined where J + = 0. Both Sui2+ and J + vanish, and 
the second term is not sufficient to offer a finite number when J + = 0, but there is 
no conceptual outrage. The last two expressions are equivalent and offer the value of 
\dz-\. 



What is 5 in J 1 oc V 5 l ? In this second order approximation, J(w) is derived from two 
relations. 

J = J_ dz_ ; dz_ = ±\/^ — (15) 

V dJ K ' 

The first equation holds since we chose the critical point z such that z Q intersects the image 

line (along the critical direction). In fact, z Q bisects the image line. The second equation 



2 A simple pendulum traces a radial arc and may not convincingly mimick a caustic curve which is a 
closed curve with interior (which can be heirachical) . If we consider a Foucault pendulum, the rotation of 
the earth should make it a closer analogy to the caustic curve as a trajectory near a cusp if not too close to 
the equator. Nonthclcss, the directional change of the velocity at a stationary point is along the radial arc 
motion, which is the motion of a simple pendulum. 
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depends on (the square root of) the ratio of Su to 8 J, and the ratio can be expressed variously 
as shown in equation (11). Thus, the title question "What is 5 in J~ l oc \/5~ l V is not a 
question well posed because the answer depends on the choice of the multiplication function. 
However, given the preferential appearance of J_ in J = dJ = J_ dz_ (since dz + = here), 
we can agree to implicitly refer to 4J_ as the multiplication function. Then, the "truly 
relevant distance" is 5 = \5u-(u )\ and is shown in Fig. 1 in comparison to the (shortest) 
distance to the caustic and "vertical distance" 3 (to the cautic curve). 

Answer: If we choose 4J_ as the multiplification factor, the distance 5 is given by \5ou2-\ 
where 5u 2 - is the normal component of 5oj = u — uj and the caustic point u a is chosen such 
that 5oj II dJ . 



-1.1 Non-preferred Caustic Point as the Origin: u c ^ uo 

If we choose an arbitrary caustic point uo c in the caustic region of interest as the origin 
of a power expansion, then 5uo = uo — uo c where u is a source position in the caustic region. 
If we let x = 5u + and y = <5cj_; u = 5z + and v = 5z_; a = J + and b = J_ for visual 
simplicity, then, the second order lens equations (61) in the neighborhood of u c are rewritten 
as follows. 

v 2 ) + b 2uv) (16) 
v 2 ) - a 2uv) (17) 
One combination of the two equations reads as follows. 

(bx - ay - 2bu) = --(a 2 + b 2 ) 2uv (18) 



x — 2u 

y 



-\ w - 

-7 ( 6 <" 2 " 



3 Given a source position oj and the caustic curve, the distance of the source position to the caustic curve 
is intrinsically well defined. However, a "vertical distance" can be defined only when a "vertical direction" 
is defined. In Fig. 6.1 of SEF and subsequently in Gaudi and Petters (2001), the "vertical direction" is 
defined by the critical direction (±E_) at a caustic point where the power series expansion is made. Since 
the choice of a caustic point as the origin of the power series expansion of the lens equation can be arbitrary, 
the "vertical direction" is not uniquely determined for a given source position u>. That is in contrast with 
the fact that the images and their magnifications are uniquely determined given a source position. In other 
words, a "vertical distance" of a source position is not an inherently meaningful quantity. In Fig. 1, the 
"vertical distance" is depicted as the distance from the source position u to the caustic curve in the critical 
direction defined at uo Q . 
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Since we have chosen a non-preferred caustic point, the image position solutions will include 
shifts along the critical curve from the origin z c where 8 J is calculated. In this second order 
approximation, the change of J± along the critical curve can not be incorporated because 
higher order derivatives of the Jacobian determinants become involved in the third order 
terms and higher. Thus, VJ is implicitly considered constant in the neighborhood of z c . 

Linear Approximation for dz + : 

We take the linear approximation as in rb99.12 (but correctly) by ignoring the RHS. 



In rb99.12, the second term (—ay) in equation (19) is missing. This missing term causes 
what amounts to an intrinsic violation of the second order nature of the equations or of the 
proper handling of a triangular matrix. Since the linear contribution for y vanishes, the 
second order contribution to y feeds back to x through the second order contributions. In 
order to cut at the linear order of 5u±, we need to fully incorporate these contributions in 
the second order in dz±. For example, if x, y, and u are an order of 10~ 3 as in a typical 
case of a main sequence source star in the Galactic bulge lensing, then v is an order of 
\/l0~ 3 and 5u2+ is of the same order as 5u>i + . In fact, we know that we can choose u a for 
the caustic point such that 5u 1+ vanishes while 5uj 2+ doesn't. Once u is determined by the 
linear equation in (19), we can find v in terms of u and y from equation (19), which is the 
same as rb99.13 with dz + reinstated. 



The first term in the RHS describes a shift along the critical curve, and so the non-vanishing 
contribution to J = d J = au + bv comes from the second term in the RHS. 



(bx — ay) = 2bu 



(19) 




(20) 




(21) 



We can rewrite it in terms of the original notations. 



J(u) 



(22) 




(23) 



If we stare at the solutions in (19) through (23) for a moment, a few things are clear. 
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l. The line (or graph) defined by (bx — ay) = 2bu in equation (19) is parallel to VJ. Thus, 
the linear approximation for dz + in equation (19) implies that the line (bx — ay) = 2bu 
should connect the source position uj and the approximate preferred caustic point. The 
approximate preferred caustic point is marked as w— in Fig. 3. The line (bx — ay) = 
2bu intersects the x-axis at x = 2u and is at a distance d from the line defined by 
(bx — ay) = as shown in the same figure. 



d 



2u 



\2dz + \ 



Va 2 + b 2 |VJ| |VJ| 



5u + — -^5lu_ 

•J 



(24) 



2. The approximate preferred caustic point determined by the line (bx — ay) = 2bu can 
be found by solving equations (20) and (30) simultaneously. 



2u 



y 



a 2 + b 2 
' -166 



a(a 2 + b 2 ) . . 9 
x « 2u - -i — -L (2uf 

m 2 v ; 

a 2 + b 2 2 a 2 + b 2 ( J + 2 

■X" ~ : U ~ : OUJ- — OW, 

Ab 166 V J_ + 



The approximations for re and accordingly for y have been made based on the validity 
condition for x in equation (29). The caustic point is marked by wO in Fig. (3) and 



will be referred to as co co in this text. 



What is Auj„ ? 

If we consider 5cu(ui co ) = uj — uj co , then Au- in equation (23) is nothing but the normal 
component of Su(u Co ). 

a 2 _|_ £2 / j \ 2 

(w - cu co ) • E„ = (5u- - y) w <5cj_ + f ^ - j 1 ^ J ( 25 ) 

We impose an implicit understanding that the caustic point cu co = u(J = 0) is to 
be found in the direction of the normal dJ calculated at z c where the power series 
expansion of the lens equation is made. 

Acj_ = (uj- u(J = 0)) • = 5uj- - 5uj^(J = 0) (26) 

If u c coincides with u C() , then u c = u co = u Q , and = 5uj_. 

3. The causic curve in the second order approximation satisfies equations in (10), and the 
resulting quadratic equation for the caustic curve is written in equation (27). If (x,y) 
is a caustic point, then x and y lie on the line defined by the first equation in (10), and 
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the line intersects the x-axis at x = 2u as shown in Fig. 3. The caustic point for the 
given value of u differs from the caustic point determined by the line (bx — ay) = 2bu, 
and the former is marked as w+ in Fig. 3. It should be clear that wO is always located 
between w+ and w—. The linear approximation for dz + in equation (19) is valid when 
w+ ~ wO ps w— with the second order approximation. 

Error from Nonlinear Correction for alz + : In order to test the goodness of the linear 
approximation for dz + , let's examine the intersection point of the second order equation 
(17) with the x-axis. If we set y — 0, equations of (17) leads to the following. 

v = J ) (~ a ^ VJ)w ; 

As we mentioned before, u is sufficiently smaller than 1 (Einstein ring radius), and the 
second term can be ignored unless it is near a cusp where b = J_ ~ or |VJ| is large 
as in the case of small caustics such as planetary caustics we will briefly discuss below. 



Quadratic Caustic Equation: Now, we write out the quadratic caustic equation and its 
validity used in the previsous section. The equation for the caustic curve can be obtained 
from equation (17) and the condition = dJ = J + dz + + J_dz- = au + bv. 

-16by= a ^-(bx + ayf =► bx = -ay ± (27) 

4. The equation is non-linear in y, but we can confirm that the origin (x = y = 0) is 
indeed a critical point [dyjdx = 0). 

2. Where by < 0, every y has two solutions. See Fig. 4. 

3. The slope dyjdx diverges where the "source line" (parallel to VJ) intersects with the 
caustic curve given by equation (17). 

y ° = - aV + & 2 ) ^ b XQO -a yoo = 0. (28) 

This indicates where equation (17) fails to describe the original caustic curve. We 
can say that the quadratic equation (17) (not necessarily quadratic in the orthogonal 
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coordinate variables x and y) is valid where \y\ « \yoo\, or equivalently where (bx) 
dominates over (—ay). From this condition, we get a criteria for \x\. 



4 




\x\ « — — - 









(29) 



For small x where |ay| << \bx\, equation (17) becomes a simple quadratic equation. 

-lQby = (a 2 + b 2 )x 2 & -16J_5cj_ = \VJ\ 2 5lo 2 + (30) 

4. If J_ = 0, the range |x| << 0. In other words, the quadratic equation for the caustic 
curve is completely invalid at the cusp. It is consistent with that the second order 
approximation of the lens equation fails where J_ ~ 0. 

If J_ = 0, the source line bx — ay = becomes y — 0, and so dx/dy = at y = 0. This 
is another way to look at the failure of the quadratic equation near a cusp. 

5. If J + = 0, the range of \x\ in equation (29) is unlimited. That is directly because 
(—ay) term vanishes and so bx term dominates. If we look at Fig. 4, bx + ay = and 
bx — ay = become degenerate, and the quadratic caustic curve in Fig. 4 becomes 
symmetric: bx = ±f(y), or —16by = (VJ) 2 x 2 . 



1. A Single Lens with Constant Shear 

In 1979, the first double quasar Q0957+561A,B (Walsh, Carswell, and Weymann 1979) 
was discovered and interpreted as two images of one quasar gravitationally lensed by a galaxy 
(or a group of galaxies) as a continuum distribution of mass. Chang and Refsdal (1979) 
considered the effect of "image splitting" of an image due to the granularity of a single 
star which happens to be near the astrometric (or transverse) position of the image. The 
angular scale of the "image splitting" of an image of a quasar (at a cosomological distance) 
by a stellar mass object is an order of microarcsecond while the two images Q0957+561A,B 
are separated by ~ 6 arcseconds. Thus, the nearity of the would-be image lensed by the 
galaxy to the (microarcsecond) lensing star amounts to the precision level of ~ 10~ 6 , and the 
gravitational influence of the lensing galaxy in the small area defined by the microarcsecond 
lensing radius of the star can be considered constant. 

The equation for a single point lens with constant shear entails a preferred direction and 
that is the direction of the (effective) mass which generates the constant shear in the small 
area of interest around the single point mass. In other words, constant shear is intrinsically 
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an approximation and can be obtained by taking large mass and large distance limits of 
more physically well-defined mass distributions. 

As a concrete example, let's take a large mass parity limit of the binary lens in Eq. (43). 
The lens equation has been normalized by Re the Einstein ring radius of the total mass, 
and we can suppose that e\ « e 2 such that e\ is the mass fraction of the (microarcsecond) 
lensing star. If we rescale the lens equation by the Einstein ring radius of the stellar mass e±, 
then ei = 1, and e 2 >> 1. We can set X\ — (choice of the coordinate origin), and then the 
microarcsecond lensed images will be around \z\ ~ 1. In order for the lensing by the larger 
mass not to dominate the lensing behavior of the images around \z\ ~ 1, the bending angle 
by the larger mass must be the order of the bending by the stellar mass €\ or smaller. For 
images around \z\ ~ 1, \x 2 \ » \z\, and the lens equation can be rewritten as follows. 

u;-^ = z--^ + e 4 + 0(x^). (31) 

X2 Z — X\ X2 

We note that ~ Re is the Einstein ring radius of the mass e 2 as a single lens. If e± is at 
a distance of ~ Re ~ from e 2 , then the source at uj ~ e 2 /x 2 produces images at whose 
positions the influence of the large mass e 2 can be considered a constant shear of order 1: 
7 = t 2 /x\ ~ O(l). 

The constant shift e 2 /x 2 indicates that the emission source is near the position of the 
large mass x 2 >> 1 for the images around with \z\ ~ 0(1), and the effect of the single 
point mass is to "split" the would-be image around the Einstein ring of the large mass. The 
outer image is splitted by a quadroid and the inner image is splitted by two trioids. These 
behavior can be very easily understood from the binary lens equation which is a physically 
well-defined closed system. 

It is customary to ignore the constant shift e 2 /x 2 in a lens equation with a constant 
shear. 

u = z — - + •yz (32) 

This equation is often referred to as Chang- Refsdal equation (e.g., SEF) and the constant 
shear term points to the positive direction of the x-axis. This particular directionality has 
been inherited from 1) the binary equation (43) where the lens axis is along the real axis 
and 2) the implicit assumption in the derviation of equation (31) that the larger mass was 
in the positive direction. In general, the constant shear coefficient 7 would be a complex 
number 7 = t 2 jx\ half whose phase angle points to the effective direction of the larger mass 
that generates the constant shear. 

It is noteworthy that the second term in equation (32) due to the (microarcsecond) 
lensing star is also a shear term, and this shear term vanishes at 00 - far away from the mass 



-14- 



that generates the shear. However, the constant shear does not vanish at oo. It is in a sense 
a reminder that constant shear is inherently an expression of approximation, and the lens 
equation with a constant shear should be treated as an incomplete equation which requires 
auxiliary assumptions or interpretations to be a physically viable model for a lensing system. 
For example, the "source position" in equation (32) can not be the position of the source 
because the "source position" is in practice the position of the would-be image generated by 
the larger mass that effects the constant shear around the (microarcsecond) lensing star. 

In order to see the effect of the constant shear, let's remove the point mass of the star 
from equation (32). 

uo = z + -iZ] J=l-| 7 | 2 (33) 

The equation is linear, and the constant shear generates only one image. The magnification 
(| J\~ r ) and parity (sign of J) of the image is the same everywhere irrelevantly of the position 
of the "source" , which is arguably the essence of a constant shear. 

z = T^f (34) 

The image is contracted in the direction of the shear 7 by a factor (1 + 7) _1 and elongated 
in the orthogonal direction by a factor (1 — / ~y)~ 1 . The distortion is universal where the 
constant shear approximation applies. In weak lensing where the shape of the objects is one 
of the measurable and interpretable quantities, such systematic distortions due to lensing 
can be detected and are being pursued to study large scale structures that may shed light 
on the history and make-up of the universe including dark stuff (matter, energy, essence, 
extra dimensions, topological defects, ... , even though it is unclear whether the essence is 
essential, extra dimensions are most likely a fundamental ingredient but for the time being 
a tantalizing yet unproven conjecture, and the exact nature of the topological defects or 
objects may turn out to be as fascinating as the existential properties of the space or the 
baryon number and remain illusive for a long time to come). How well one can determine 
the shear distribution crucially depends on the shape resolution and statistics of the shape 
objects. A degree scale large format space imager Galactic Exoplanet Survey Telescope will 
be a fitting explorer even though the wings have been clipped again. 



1.1. The k Field of the Single Lens with Constant Shear 

The linear differential behavior is insensitive to the constant shift of the source position 
and can be discussed routinely as usual. As in the case of n-point lenses (Rhie 1997; Rhie 
& Bennett 1999; Rhie 2001), the linear differential behavior of the lens equation (31) 
is completely determined by one analytic function k, and so we use and follow the facts, 
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conventions, and analysis patterns in the references (RBCPLX from here on). The Jacobian 
of the lens equation can be written as follows where the 2x2 matrix indicates that the 
underlying space of the complex plane is a 2-d linear space. 

J=( Kl ); « s «s s ^ = F+7 < 35 > 

The condition for criticality, \k\ — 1, is a second order analytic polynomial equation, and so 
the phase angle of k changes by 2n along the critical curve. 

k = e i2lp ; Ay? = 2tt (36) 

Thus, the total topological charge of the critical loops or caustic loops is ^2\e\ = 1, and 
the caustic is made of one 4-cusped quadroid or two 3-cusped trioids. The depiction can be 
found in Chang (1984) and also in SEF.Fig.8.8. There are two limit points (where k — 0) 
on the imaginary axis (orthogonal to the direction of the constant shear) at z — ±2/^/7, and 
the trioids enclose the limit points. 

J(z) = 1 - \k(z)\ 2 < 1 ; J(oo) = l-7 2 (37) 

The Jacobian determinant J(z) of this lens equation with a constant shear is bounded by 
J(z) < 1. The maximum holds at the finite limit points z = ±2/^/7 but not at infinity. In 
fact, the critical curve can pass through the infinity and the image at infinity can be infinitely 
bright. This unphysical behavior derives from that the lens equation with a constant shear 
is accommodated by an infinitely large mass at infinity when the lens equation is extended 
to the whole lens plane. This behavior is easy to understand when seen as an approximation 
of the binary lens. The critical points at infinity correspond to the bifurcation points where 
the topology of the critical curve changes from one loop to two loops. At the same time, the 
topology of the caustic curve changes from one quadroid to two trioids. 



1.2. Images 

We have derived the lens equation (32) of a single lens with a constant shear from the 
binary lens equation and can read off some of the basic properties of the equation following 
the recipes developed for binary (or n-point) lenses briefly reviewed in the following section. 

1. The caustic loops are simple loops. In other words, they neither self-intersect nor nest. 
Thus, the source plane accomodates only two domains: V° (outside) and V 1 (inside). 

2. There are two images outside the caustic loop(s). One only needs to examine the 
number of images of uj = 00 because the outside domain V° includes 00. 
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3. A source inside a caustic loop (V 1 ) produces four images, and the caustic loop itself 
produces three image loops one of whose is the critical curve. The caustic curves have 
cusps, and so the two "non-critical curves" have kinks at the corresponding points 
unless they are precusps. However, the critical curve is smooth because the inverse 
mapping of the lens equation from the caustic curve to the critical curve is singular at 
the cusps the singular points. Precusps are bifurcation points of the critical curve and 
"non-critical curves" (see Fig. 10 in Rhie (2001)). 

The constant shear breaks the axial symmetry of the single lens but leaves a residual 
reflection symmetry as is familiar from binary lenses (Rhie & Bennett 2001). If the source 
is on the lens axis, u = uj, there are two images on the the lens axis, z — z, and two extra 
images on a circle l^ -1 = ^/l — 7 : 7 > 1. The lens equation subject to uj = Q leads to the 
following equation. 

= (^)(l-^-7) (38) 

1. z = z: The two images on the lens axis are very much like the two images of a single 
lens that form on the line defined by the positions of the lens and source. 

u = z{l+ 1 )-- =>■ U = ^1 + 72 - - (39) 



z 



[\/r+7 z\ 



The second equation shows the correspondence with the equation of a single lens. 

2. \z\~ 2 = 1 — 7 (> 0): The two images are on the circle of radius \z\ = l/y/1 — 7 
centered at the lens position (here the origin) where the real component of the image 
positions is given by a simple linear scaling of the source position from the lens. 

u = 2jRe(z) : 7 > 1 and \Re(z)\ < —=L (40) 

V 1 - 7 

When 7 < 1, the caustic curve is a quadroid intersecting the lens axis. The cusps of 
the quadroid are easily calculated: two are on the real axis (the lens axis) and has 
k — 1, and the other two are on the imaginary axis and has k = — 1. 

±2 7 ±i 2 7 

The two images on the circle exist when u is inside the quadroid, which can be easily 
read off from equations (40) and (41). The Jacobian determinant of the images on the 
circle is given as follows. 

ie 

J = 47(1 - 7) sin 2 6 >0 : z = - (42) 

V 1 - 7 
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The images on the circle are positive images except at the cusps. At the cusps on the 
real axis given in equation (41), J = 0, and the circle is tangent to the critical curve. 
We can easily check that the images on the real axis are inside the critical curve and 
so are negative images. The total parity of the images of a single lens with a constant 
shear is zero. 

When 7 > 1, the caustic curve is made of two triods enclosing the two limit points, 
and the two caustic loops are off the lens axis. Thus, u = u has only two images and 
they have opposite parities exactly as in the case of a single lens. 

The fact that the two extra images are on a circle deserves a modest attention. If we 
decrease the effect of the constant shear, 7 — » 0, the circle approaches the critical curve of 
the single lens of radius 1, namely the Einstein ring radius: J — > and \z\ — > 1. In general, 
Einstein ring, critical curve, and image ring are all distinct. Einstein ring determines the 
(transverse) distance scale of lensing. Critical curve is where the Jacobian determinant of 
the lens equation vanishes, images are degenerate, and the magnitude of the images diverge. 
Image ring is a proxy ring which a finite size source fallen inside a caustic loop produces, 
and its shape and intensity depends on the caustic structure and the relative position of the 
source. In the case of a single lens, they all converge to one ring. In the case of the single 
lens with a constant shear, its proximity to a single lens makes the image ring quite circular. 
In fact, the same line of thought offers a valid intuition for caustics of higher multiple lenses 
when the effect of the lens elements but one can be considered a long distance effect. 



2. Binary Equation, Quasi-Analytic Equations 

The binary lens equation is written as follows (see RBCPLX, references therein, or 
others). 

to = z ^— , (43) 

Z — X\ z — x 2 

where u, z, and Xj : j — 1, 2 are the positions of a source, an image, and the lens elements 
of fractional masses e 3 ■, : j — 1, 2 in the two dimensional sky as a complex plane. We have 
chosen the lens axis to be along the real axis so that x\ and x 2 are real, and the unit distance 
scale is given by the Einstein ring radius Re- 

1 — Re — V&GMD ; ^ = ^ + ^ (44) 

where D is the reduced distance. 

When there are n point lens elements where n > 2, we only need to extend the range 
of the index to j = 1, 2, n and replace the lens position vectors to Xj because we can not 
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line them up all on the real axis. The form of the Jacobian matrix in equation (35) applies 
to any quasi-analytic lens equation, and k = dcu is easily calculated for each class of lens 
equations. 

K = V" e ± 



J 



(z - Xj y 



J = 1 - \ K \ 2 < 1 ; J(qo) 



2.1. Eigenvalues, Eigenvectors, and Orthogonal Decompositions 

The quasi-analytic lens equations - the lens equation is not analytic, but the differentials 
are - have been discussed in RBCPLX and is reviewed here briefly. The eigenvalues of the 
Jacobian matrix are easy to find, and A_ vanishes on the critical curve (\k\ — 1 44> J = 0). 

A± = 1 ± \k\ , (45) 

If 2(p is the phase angle of k such that 

k = \n\e 2iv , (46) 
the eigenvectors are given by (±)e± where we choose the basis vectors e± as follows. 

V2 , (47) 



(+ = l&J' e -=\-ie^) ; l|f=l 

The orientation of the basis vectors defined here can jump or mismatch around the 
critical curve. Here we are mostly concerned with the local behavior of the images near the 
critical curve, and the global orientability of the basis vectors as chosen in (47) along the 
closed curves would not be a relevant concern. 

The complex plane is a real plane with complex structure, and it is sufficient to write 
out half the 2-d equation because the other half is the complex conjugate of the former half. 

duo = dz + Rdz + \dMz 2 + 0{\dz\ 3 ) (48) 

We only need to consider the upper components of the basis vectors, which are mutually 
orthonormal complex numbers. 

E+ = e~ iLp ■ E_ = ie~ llp (49) 

(E + ,E-) forms a right-handed coordinate frame, and any vector can be decomposed into 
the orthogonal components. For example, dz = dz + E + + dz_E^ where the linear coefficients 
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dz+ and dz^ are given by the inner product {dz, E±) that is uniquely defined from the norm 
of a complex number. If dz = \dz\e ld and r = \dz\, 

dz + = Re(dzE + ) = r cos(# + <p) = rC (50) 
dz_ = Re(dzE_) = r sm(9 + (p)=rS (51) 

We note that the decomposition coefficients are real. 

We recall from RBCPLX that E + is always tangent to the caustic curve and so E_ is 
always orthogonal to the cautic curve. Thus, decomposing the quantities relevant in lensing 
into orthogonal components in E + and is a bit more than an idle exercise. The Jacobian 
determinant of the lens equation J(z) = 1 — |/t| 2 is a real function and its differential dJ can 
be calculated in any basis we are pleased to choose. If the 2-d real plane is parameterized 
by (£1,^2) such that z — £1 + i£,2 and z — £1 — ^2, then d = d z — (d\ — 182) /2 and 
d = dz — (di + id 2 )/2. We may call these coordinate changes complexification and express 
z = (£i,6) C and 2d=(d u d 2 ) c . 

dJ = d 1 Jd£ 1 + d 2 d& 
= dJdz + dJdz 
= d + Jdz + + d-Jdz^ 

The gradient VJ = (diJ,d 2 J) points in the direction of the maximum increase of J and is 
a normal vector to the critical curve. So is dJ. 

(VJ) C = 2dJ (52) 

From equation (52) and conversion between (dz + ,dz_) and (dz,dz), we can confirm that 
d + J and <9_ J are -E± components of the normal vector (VJ) C . 

2<9J = d + JE + + d„JE„ (53) 
= J + E + + J_E_ (54) 

We remind that J + and J_ are real, and so are d + and cL. 

9± = ^±9 + E±d (55) 

From here on (and retroactivley if applicable) we freely interchange between VJ and 2dJ 
because they are one and the same vector expressed in different coordinates. 
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2.2. Near the Critical Curve and Caustic Curve 



The variation equation (48) can be sorted out into E± components. In the linear order 
in r = \dz\, 

du> = rC(l + \K\)E + + rS(l- \k\)E_ , (56) 

and the i?_-component vanishes on the critical curve because \k\ = 1. This vanishing linear 
component is behind the quadratic behavior of the lens equation in the critical direction 
(E-) in the neighborhood of the cirtical curve and is at the foundation of the square-root- 
of-the -distance dependence of the magnification of the images near the critical curve. We 
have discussed fully in section -1 that the lens equation becomes a real symmetric quadratic 
equation from the image line (±i?-) to the source line (3 J) where J /a 0. 

We calculate the second order term in equation (48) at the critical point. 

\dn dz 2 = --r 2 ((J+C 2 + JS 2 )E + + (J_C 2 - J+S 2 )E.) , (57) 

where C 2 = cos 2(0 + </?) and S 2 = sin 2(6* + </?). The E ± -components of the linear deviations 
and second order deviations can be listed as follows. 

5u 1+ = 2rC (58) 
5cji_ = (59) 

8u 2+ = -^r 2 (J + C 2 + JS 2 ) (60) 

<^ 2 _ = -I r 2 (J_C 2 - J + S 2 ) (61) 

If we consider a unit vector unit = (C 2 , S 2 ), then 5u 2+ depends on the unit-component of 
VJ (inner product, or dot product), and 8u 2 - depends on the signed area (rotation, or cross 
product) defined by unit and VJ. It is usually the case 4 that VJ 7^ 0. The third order 
terms have the similar structure as the second order terms. 

5u 3+ = -lr 3 ((J ++ - J__)C 3 + 2J + S 3 ) (62) 

^3- = -^r 3 (2J + _C 3 -(J ++ - J-)S 3 ) (63) 

The reason is because K-field is analytic and so is orientable. In the third order, the phase 
is multiplied three times, and so the unit vector is replaced by unit = (C 3 , S 3 ). 



4 The gradient vanishes on the critical curve only where the critical curve bifurcates and changes its 
topology. In the case of a single lens with constant shear 7, there are two bifurcation points at 00 when 
I7I = 1. If we go back to the derivation from the binary lens equation, |— y| = 1 holds when the constant shear 
is calculated on the Einstein ring of the larger mass. 
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2.3. Rough Estimation of |VJ| 

In a typical Galactic bulge lensing where a star in (the direciton of) the bulge is lensed 
by a (faint) star or a planet system in the bulge or the disk, the stellar radius in units of the 
Einstein ring radius of the total mass of the lensing system is r* ~ 10 -3 — 10~ 2 . 

In order to get an idea of the magnitude of the gradient J on critical curves, let's consider 
the simplest case of a single lens. For a single lens, k = 1/z 2 ; J = 1 — 1/|-2| 4 ; |VJ| = 4/|z| 5 . 
On the critical curve, \z\ = 1, and |VJ| = 4 >> r*. The gradient J is normal to the ring 
\z\ = 1 and so is the positive eigenvector (±)E + . Thus, J_ = 0, and J + = (±)4. The fact 
that J_ = everywhere on the critical curve underlies that the point caustic is a degenerate 
cusp. 

In a binary lens where its lens elements masses are comparable, we expect a similar 
range of values of |VJ|. We have shown the numerical values of a/|J-| on the 4-cusped 
central caustic of a binary lens with the mass ratio of ~ 2 : 1 in Fig. 1 of Rhie & Bennett 
(1999) (rb99). We also stated that trioids have somewhat larger values. The larger the value 
of |VJ| around the critical curve, the faster decreases the magnification value and smaller 
the "width" of the caustic curve. 

If the mass ratio is very small as in a planetary binary lens, then 

(V J) c = 2d J = -2kBk = 4kV - €j _ Nq , (64) 

i ( z ~ x j) 



and on the critical curve of a small mass lens e 2 « 1, 

IVJI R 



£ 2 



4 , x 

>> 1 . (65) 



For terrestrial planets, the mass ratio to the host star is ~ 10~ 5 — 10~ 6 , and the second or 
high order terms can be dominant whether the planet is bound or free-floating. Thus, power 
expansion approximation can not be used for integrating over the stellar luminosity profile. 
In a multiple point lens system, the range of | VJ| can be very large depending on the caustic 
loops. It is safe to say that the central caustic will have | VJ| ~ 0(1). In the practice of light 
curve fitting, model light curves are all calculated numerically. At caustic crossings, what 
may be referred to as "local ray-shooting method" is used to be able to accomodate the rapid 
change of the inverse Jacobian determinant J' 1 and instability related to the criticality. 
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-2. Conclusion and Discussion 

We have reexamined the quadratic behavior of the lens equation near the critical curve. 
In the second order approximation, the lens equation becomes a real symmetric quadratic 
equation between the source line and image line locally defined at a proper causic point and 
critical point respectively. In other words, given a source position uj, there exits a preferred 
caustic point uj such that 5uj = uj — uj is in the same direction as dJ(z ) where z Q \— > uj 
under the lens equation. The source uj generates two images at z± along the critical direction 
defined at z , and the magnification of the images is given by a/ (45oj_ J_) _1 . Source line 
refers to uj — uj q and image line refers to the line connecting z± bisected by z . We emphasized 
that the image line (±E_(z )) and source line (dJ(uj )) are not orthogonal to each other. 
In fact, the second order approximation fails when they are orthogonal because they are 
orthogonal at cusps (J_ = V J • E_ = 0). 

It is not necessarily a trivial matter to find the preferred caustic point uj for each 
source position uj, say, of a time sequence of a stellar disk that crosses the caustic curve. It 
is in fact unncessary because we can recover the image positions and their magnifications by 
power expanding the lens equation at a nearby but arbitrary caustic point. Fig. 1 clearly 
indicates that we have freedom to choose a caustic point which is in the neighborhood of the 
caustic point that defines the (shortest) distance. 

The criticality (vanishing derivative) occurs only in one direction, and the power ex- 
pansion of the lens equation at a non-preferred caustic point (meaning ^ uj Q ) involves non- 
vanishing non-critical component of dz (i.e., dz + ^ 0) and this non-critical component is 
a linear equation of the source position variable components uj + and 5uj_ in the lowest ap- 
proximation. The non-critical component dz + causes a migration along the critical curve 
and we can see in Fig. 3 that the non-critical component 2dz + measures the shift of the non- 
preferred caustic point from the preferred caustic point in the direction of E + . If we include 
the nonlinear effect, the shift deviates from 2dz + . 

In this Taylor expansion approximation, the gradient of the Jacobian determinant 
(VJ = 2dJ with the understanding that we freely exchange between complex notations 
and real countparts for the same 2-d vectors) is a constant vector calculated at the critical 
point where the power expansion is made (or at the origin). And, the directional difference 
of d J at two different caustic points (uj c and uj cq ) can not be incorporated in the second 
order approximation. In the lowest approximation, dz + is a linear equation of 5uj + and uj_, 
and so the direction of the source line of uj is given by the source line of a source position 
whose preferred caustic point would be the origin (uj c ). The preferred caustic point for uj 
in this approximation is given by what we denoted uj co ) (and marked as wO in Fig. 3). We 
emphasize that the preferred caustic point uj co ) is distinguishable from what we referred to 
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as the approximate preferred caustic point that is determined as the cross section of the 
linear equation bx — ay = abu and the quadratic caustic equation. The linear approximation 
for dz + is valid Where the difference is sufficiently small. 

Once we know the preferred caustic point, the gradient of the Jacobian determinant, 
and the normal to the caustic where both of the latter two are calcualted at the origin (a> c ), 
then the image positions and their magnifications can be calculated. The magnifications 
formula can be written in the same form as those calculated at the preferred caustic point 

J = ±v/J_Acj_ (66) 

Ac<j_ is the normal component of the shift of uo from the preferred caustic point cu CQ . We 
illustrated in Fig. 3 that normal component and "vertical distance" 5 are clearly distinguish- 
able both conceptually and practically, hen the curvature of the caustic can be neglected 
either because the stellar disk is sufficiently small or J + ps in the neighborhood of the 
caustic crossing point, both the "vertical distance" and the (shortest) distance become com- 
parable to the normal component. When u CQ and u c coincide with oj , Acj_ = 5u-, and this 
symmetric quadratic case is depicted in Fig. 1. 

J = ± v / XaZ (67) 



Gaudi & Petters (2001) raised an issue of the identity of the "relevant distance". As 
we have emphasized many times, the "truly relevant distance" is |Au;_| given as follows 
(rb99.16 but with the proper substitution of the linear approximation of dz + in terms of 5u + 
and 5u)-). 

\dJ\ 2 

Au- = 5uj- - 8u-(J = 0) = 5u- + —r^dz% , (68) 

U 

where the linear function is given as follows. 

(J^5uj + - J + 5uj_) = 2J_dz+ (69) 

We reiterate that 5uj_ is the E_ component of 5u = uj — uj c where uj c is the origin, and 
5uj_(J = 0) is the E_ component of the position vector of the caustic point (ou co ) which lies 
on the source line. The criticality dictates two normals, and the lens equation in the second 
order approximation is a quadratic relation between the variables along these two normals. 
And, the two normals are oblique to each other. 



5 According to The Random House College Dictionary, vertical is to be in the same direction as the axis. 
Here the axis is taken as the critical direction vector at the origin (E- (u) c )) where the power series expansion 
is made. The critical vector at the origin is normal to the caustic curve at the origin. 
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It should be worth pondering if the notation 8uj-{J = 0) is potentially confusing. What 
else can one imagine for the quantity 8iv( J = 0) aside from the intrinsically relevant quantity 
that is the derived preferred caustic point position vector uj cq 1 One can drop a straight line in 
the direction of the "vertical axis" from uj in Fig. 3 and claim that the intersection point with 
the caustic curve must be 8uo(J = 0) and its E_ component be 8u-(J = 0). That may seem 
all right at a first glance. But, what would be the significance of the line that is parallel to 
E_{u) c ) and passes through ul We do not know. However, we know that the corresponding 
equation 8u + = constant (given by the value of the E + component of the particular source 
position) can not be derived from the quadratic lens equation (17). Given the irrelevance, we 
can discard the case from the list of potentially confusing interpretations. What else can one 
imagine for 8u(J = 0)? Currently, we lack imagination for other possibilities of confusion. 
We take it as a good enough reason to pardon our notation and close the case. That is, of 
course, until someone brings a brilliant confusion candidate to our attention. 



1. The interpretation of Gaudi & Petters (2001) that |Ac<j_| is the "vertical distance" is 
an erroneous misinterpretation. 

2. Gaudi & Petters (2001) claims to have reproduced the results in SEF. We find that the 
second equation of SEF. 6. 17 has a mistake. By now, we have a clear understanding that 
when we expand the lens equation in the neighborhood of an arbitrary non-preferred 
caustic point, there will be terms describing a migration along the critical curve. The 
condition for migration is = d J = au + bv, and so we expect v = —au/b + .... If 
the two equations in SEF. 6. 17 are compared (x\ is correct), it should be clear that the 
missing term is a 2 y/(2b 2 ) in our notation. 

(a) The fault is at the truncation prescription in the first paragraph in SEF. p. 189. 
If we look at the RHS of equation (12), the most dominant term is v 2 term, and 
so v ~ yfy. We know that that is exactly the case, namely v oc y/y, when the 
origin is the preferred caustic point and so there is no migration along the critical 
curve. If we look at equation (11), it is clear that the dominant term in the RHS 
is v 2 term and so u is linear in x and y as is shown in equation (14). Thus, 
v 2 ~ x ~ y ~ u, and they are sufficiently smaller than 1 (Einstein ring radius). 
Incidentally, that is why v is much larger than the other three, and images across 
the critical curve are tangentially stretched large. Since the migration condition 
is linear in u and v, we expect v to be made of terms of two different orders: ^fy 
and linear terms in x and y. Therefore, all the terms up to the second order in y 
(so, xv ~ v 3 ,xv 2 ~ v A ) should be kept to calcualte v from the second equation of 
SEF. 6. 17. The linear terms in x and y come from xv and v 3 terms respectively. 
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(b) One may suggest to ignore the linear term in y because that is smaller than yjy 
term. One problem is that the linear term in y does not have any reason to be 
small in comparison to the linear term in x, which is reflected in the first equation 
in SEF.6.17. In other words, there is no reason to truncate v 3 term while keeping 
xv because they are of the same order. Thus, there is no foundation for the 
prescription to ignore the linear term in y while keeping the linear term in x. In 
fact, dropping the linear term in y causes a conceptual inconsistency by obscuring 
the geometric nature of the linear terms - migration along the critical curve. One 
may argue that a = J + ~ 0, of course. Then Au_ approches 5u_, and the 
(shortest) distance and "vertical distance" converge to the normal component. In 
general, we should keep the both linear terms. 

The missing term a 2 y/ {2b 2 ) affects the expression of Ac<j_. If one cosmetically interprets 
the erroneous expression, one can be led to the erroneous conclusion on the distance 
drawn by Gaudi & Petters (2001). The fault may lie in that the apparent relevance of 
the vertical distance which may have been a pleasing discovery for a brief moment has 
not been tested through the routine logical digestion process necessary for a meaningful 
claim and has been put forward as a fact. The correct relevant distance is the normal 
component 5uj^ — 5uo_(J = 0) where the normal direction is given by the critical 
direction at the origin where the functional definition of the origin derives from that 
that is where the Taylor (or power) expansion is made. 

3. The error in rb99 we described somewhat extensively in a previous section is propagated 
to rb99.13 through rb99.16. The substitution 5u + <= 2dz + in those equations should 
be replaced by the correct substitution Su + — 5u-J + /J- <= 2dz + . 



-3. Epilogue: Another Long Walk to Renaissance 

Gaudi & Petters (2001) emphasized to have worked out in (2-d) real variables. They 
made an erroneous conclusion and kindly warned the readers accordingly, 'there are not 
many people who know that the distance is the "vertical distance" . ' This claim may imply a 
correlation between the error and the choice of variables, but that is not the case. The missing 
term in the second equation of SEF.6.17 could have been recovered if one had stared at the 
two equations of SEF.6.17 and scratched the head about the meanings of the linear terms 
in comparison to the square root terms. In our humble belief, (astro)physically meaningful 
quantities have intrinsic values (many times as the form of geometric quantities) and offer 
intuitive understanding. 
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One advantage of the complex variable is that quasi-analytic lens equations are, well, 
quasi-analytic. The derivatives are completely determined by one analytic function. There 
is nothing much great about it except that it is simpler. The system is constrained and the 
constraint is completely contained within the nature of the variable, and so it is easier to 
get to know the system. Especially, for one variable analytic function, we even know how to 
integrate because the first theorem of analysis is simply defined. 

Is there a taboo against complex variables in astro-physics? If there were, it must be 
the time to change now that astro-physicists are gearing up to the idea of testing red-shift 
dependence of the electromagnetic fine structure constant - the beloved prime number 137 
when inversed and truncated. (S. Bechwith, HSL workshop, U.Chicago, April, 2002; See 
J.K.Webb et al, PRL, 87, 091301 (2001) = astro-ph/0101375 for a ~ 2a level claim of 
Aa/a ~ — 10 -5 from quasar absorption spectrum analyses 6 ; many more papers in the usual 
suspect sites astro-ph and hep-ph; P.A.M.Dirac may have been the first to be serious enough 
to write a paper, Nature, 139, 323 (1937), on the variability of the fundamental constants 
including the fine structure constant; the first Kaluza-Klein models were written in the mid- 
1930's as well; see T. Chiba, gr-qc/0110118 for a review of constraints on the variation of 
fundamental constants; coupling constants are believed to run with energy; there are strong 
experimental indications for grand unification of the electroweak and strong interactions (see 
F. Wilczek, hep-ph/0101187 for a summary) and relevance of supersymmetry 7 .) 

Gravitational lensing is an old subject revitalized with sensitive and gigantic instru- 
ments. Leonardo da Vinci seems to stand out as a Renaissance man. And, Mona Lisa a 
Renaissance woman. Mona Lisa is believed to be da Vinci himself sane frizzy hair and beard. 
The mysterious smile of Mona Lisa may be a smirk of da Vinci daring the admirers, "You 
are looking at me." Or, "Can we live without phases?" It is the time to release complex 
variables from the hit list and allow some freedom. In fact, eigenvectors, eigenvalues, and 
all that are about the differential behavior in the underlying 2-d space. E± decomposition 
coefficients are all real. It is inevitable to deal with two real variables in one way or another. 



6 We do not pretend to have understood the systematics related to the analyses. Some useful numbers 
can be found in Table 1 of M.T. Murphy et al., astro-ph/0012421 and they indicate that the usage of the 
rest frame UV lines requires a precision level of <~ 10 ~ 8 in the spectral measurements to be able to test 
Aa/a ~ — 10 -5 . The resolution of ~ 7km/ sec obtained from HIRES/Keck I is far worse than the necessary 
precision, and the analyses rely on broad spectral profile fitting and x 2 estimation. 

7 We happen to be a believer of supersymmetry because that seems to uniquely naturally solve the "missing 
spin 3/2 problem". No fundamental spin 3/2 particles have been seen yet. Of course, Higgs particles have 
not been found either. Once named God particles, they may be shy away from the physical nature. LHC 
will tell us within its limit, no doubt. 
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Then, what significance does it carry to claim to have chosen real basis vectors (implitly as 
two column vectors (1,0) and (0,1)) instead of E±l It is unclear. The authors Gaudi & 
Petters (2001) did not specify the qualification of the statement for one purpose or another. 
It remains to be a mysterious invitation to a dark age. 
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Fig. 1. — Quadratic behavior of the lens equation near a critical point z Q where J-(z ) 56 0: 
The criticality defines the image line dz \\ ±E^(z ) and the source line 5u \\ dJ(z Q ), and 
the lens equation becomes a real symmetric quadratic equation from the image line to the 
source line: 5u = \d.J dz 2 _. There are two images when 5u/dJ > 0, and no images when 
5u/dJ < 0. Here 5u || dJ, and u> generates two images z± near the critical point z . The 
positive image at z + is in the direction of V J, and the negative image at z_ is in the opposite 
direction. The amplification of the images is a/ (45u;_ J_) _1 , and both the shortest distance 
and the "vertical distance" differ from the "truely relevant distance" \8u-(u )\. 
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Fig. 2. — If we consider a large stellar disk, the source line (8 J) for each point in the disk 
can be diverse. Given a source position to inside the caustic curve, u will define the source 
line, and the image line of the two images will be at z Q : z Q i— > u) Q as we have discussed in 
the text and the caption of Fig.l. One issue is how well we can recover the images z± and 
their amplification ■sj (45 ui- J-) -1 when we take a nearby but arbitrary caustic point uj c as 
the origin where the coordinate system {E + ,E_} is defined and the power series expansion 
coefficients are calculated. 
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Fig. 3. — Non-preferred caustic point oj c as the origin. The value of multiplication function 
J_ is calculated at the origin, and Au;_ is the normal component of Sou — 5u(J = 0). The 
bottom panel compares the caustic points: w+ is the correct caustic point for given 2dz + ; w— 
is the caustic point defined by the line bx — ay = 2bu and was referred to as the approximate 
preferred caustic point for u; wO is the derived preferred caustic point for u in the linear 
approximation of dz + . 
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Fig. 4. — The caustic curve in the second order approximation satisfies a quadratic equation 
between y and bx + ay. Here f(y) = a/— lQyb 3 / (VJ) 2 . The equation is a valid description 
of the original caustic curve in the neighborhood (here circled for illustration) where \ay\ 
is negligible in comparison to the other two terms bx and f(y) and becomes a quadratic 
function between y and x. 



